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The Chalker Coddington quantum network percolation model is numerically 
pertinent to the understanding of the delocalization transition of the quantum 
Oh Hall effect. We study the model restricted to a cylinder of perimeter 2M. 

We prove firstly that the Lyapunov exponents are simple and in particular that 
the localization length is finite; secondly that this implies spectral localization. 
Thirdly we prove a Thouless formula and compute the mean Lyapunov exponent 
'— 1 which is independent of M. 

1 Introduction 

^ We start with a mathematical then a physical description of the model. Fix the 

parameters 

O r,t G [0,1], such that, r 2 +t 2 = 1, 

O 



denote by T the complex numbers of modulus 1 and for q = ((71,(72, 93) G T 3 , 
by S(q) the general unitary U(2) matrix depending on these three phases 
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Let (n,T,P) be the probability space: := (T 6 ) (2Z)2 , F := <g) (22) 2d 6 Z where 



dl is the normalized Lebesgue measure on T, and F the a— algebra generated 
by the cylinder sets. With 

p G p(2j, 2k) =: ( pi,P2,P3 , P4,P5,P6 ) 

p e (2j,2fc) Po (2j+l,2fc+l) 

and the basis vectors e M (p) := (p, p G Z 2 ), the family of unitary operators 

U{p) : 1 2 {I?) -> Z 2 (Z 2 ) 
is defined by its matrix elements U^ v = {e^JJe v ): 

U^-u '■= except for the blocks 



[ ^(P)(2j+l,2fc);(2j,2fc) ^(P)(2j+l,2fc);(2j+l,2fc+l) 

\ ^(P)(2j,2fc+l);(2j,2fc) ^(P)(2j,2fc+l);(2j+l,2fc+l) 

^(2j+2,2fe+2);(2j+2,2fe+l) ^(2j+2 ) 2fc+2);(2j'+l,2fc+2) 

^(2j+l,2fe+l);(2j+2,2fe+l) ^(2j+l ) 2fc+l);(2j+l,2fc+2) 



S( Pe (2j,2k)) (1) 



S(p (2j+l,2fc+l)). 



Note that U is an ergodic family of random unitary operators; indeed, 
U*U = I = UU* because of the unitarity of the blocks; further denote by 
G the action of 1? on functions / on 1?: 

(%m)/)0) := /(A* + (2Z, 2m)) (/i G Z 2 , (Z, m) G Z 2 ), 

and, by abuse of notation, the corresponding shift on Q. Then is measure 
preserving and ergodic on SI and 



U(Qp) = QU{p)Q- x . 

This model was introduced in the physics literature by Chalker and Cod- 
dingtion, [CC], see [KOK] for a review, in order to study essential features of 
the quantum Hall transition in a quantitative way. U describes the dynamics of 
a 2D electron in a strong perpendicular magnetic field and a smooth bounded 
random electric potential which is supposed to have some array of hyperbolic 
fixed points forming the nodes of a graph. 

In this picture the electron moves on the directed edges of the graph whose 
nodes are "even": {(1/2,1/2) + (2j,2k),j,k G Z} or "odd": {(1/2,1/2) + 
(2j + 1,2k + 1), j, fc G Z} with edges connecting the even (odd) nodes to 
there nearest odd (even) neighbors. U describes the evolution at time one of 
the electron. The edges are labeled by their midpoints. They are directed in 
such a way that U models the tunneling near the hyperbolic fixed points of the 
potential, see figure 1. The tunneling is described by the scattering matrices 
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Figure 1: The network model with its incoming (solid arrows) and outgoing links 

S associated with the even, respectively odd, nodes. The i.i.d. random phases 
associated with each node take into account the deviation of the random electric 
potential from periodicity. 

Following the literature on tunneling near a hamiltonian saddle point, [FH], 
[CdVP], the parameter t is v /jq^g- where e is the distance of the electrons energy 
to the nearest Landau Level. An application of a finite size scaling method to 
their numerical observations led Chalker and Coddington [CC], see also [KOK], 
to conjecture that the localization length diverges as t/r — > 1 as 




where the critical exponent a exceeds substantially the exponent expected when 
a classical percolation model is applied to the problem, [T]; the values advocated 
for a are 2.5 ± 0.5 for the quantum and 4/3 for the classical case. 

Because of its importance for the understanding of the integer quantum 
Hall effect the one electron magnetic random model in two dimensions was and 
continues to be heavily studied in the mathematical literature. Mathematical 
results concerning the full Schrodinger Hamiltonian can be traced from the 
following contributions and their references: [W] for percolation, [GKS] for the 
existence of the localization-delocalization transition [ASS, BESB, G] for the 
general theory of the quantum Hall effect. For results concerning a 2D electron 
in a magnetic field and periodic potential, which corresponds to the absence 
of phases here, see [TKN 2 ], [HS]. For recent work on Lyapunov exponents on 
hamiltonian strip models see [RS], [BS], [Bou]. 

Our results concern the restriction of the model to a strip of width 2M and 
periodic boundary conditions; they are presented as follows. In section 2 we 
analyze the extreme cases, r = and r = 1. Then, for the case where all 
phases are chosen to be 1, we give a description of the spectrum. Questions 
related to transfer matrix formalism are handled in sections 3, 4, 5 . In section 6 
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we prove simplicity of the Lyapunov spectrum and finiteness of the localization 
length. In section 7 we prove a Thouless formula and show that the density 
of states is flat which implies our results on the mean Lyapunov exponent. In 
section 8 we prove complete spectral localization. 

2 Some properties of the model 

2.1 Extreme cases 

Note that in case of complete "reflection" or "transmission" the system localizes 
completely: 

Proposition 2.1 Let rt = 0. Then, for any p £ Cl, the spectrum ofU(p) is 
pure point. 

Proof: Assume r = 0, p G Cl and define the family of subspaces {T~ij^)(j,k)ei 2 
as: 

Hj,k = Ra n (e2j,2fc, e2j+l,2fcj ^2j+l,2k-l, ^2j,2k-l) ■ 

These subspaces are invariant under U(p) and 

®(j,fc)ez 2 Hj,k = l 2 {^ 2 ) , (2) 

which means the operator U(p) is pure point. The case t = is treated similarly. 
■ 

On the other hand one has complete propagation if all the phases are equal 
to one; define Cl 9 p = (. . . , 1, 1, 1, . . .) by p(2j, 2k) := (1, 1, 1, 1, 1, 1) then we 
have: 

Proposition 2.2 Letrt ^ 0. Then, the spectrum ofU(. . . , 1, 1, 1, . . .) is purely 
absolutely continuous. 

Proof: We make use of a decomposition similar to (2) and define the unitary 
V from l 2 (l?) to 1 2 {I?) C 4 by Ve 2jt 2k ■= e jtk (8) e lt Ve 2j+h2 k+i ■= &j,k ® 
e2,Ve 2 j,2k+i ■= ej,k ® e3,Ve 2 j + i j2 k := e.j,k ® ^4- Let P be the projection 
P := I <8> (|ei)(ei| + | e2) {e2 1 ) ■ From the definition of U in (1) one reads that 
VU 2 V~ 1 commutes with P and that PVU 2 V~ 1 P is equivalent to 

/ rt(T 0>1 - Ti, ) r 2 T lfi + t 2 T 0A \ 
V t 2 T ,_i + r 2 T_ li0 ri(T_ 1)0 - T _i) J 

with the translations on / 2 (Z 2 ) defined by 

T n ,mip(j, k) := ip(j + n,k + m) (n,m£Z). 
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The Fourier transform F : l 2 {1?) — > L 2 (T 2 ) transforms the translations to 
multiplication operators: TT n ^ m T~ x = exp (— i(nx + my)), thus the restric- 
tion to the range of P of FPVU 2 V~ 1 PF~ 1 is equivalent to a matrix valued 
multiplication operator 

/ ri(e-^ - e~ ix ) r 2 e~ ix + t 2 e~^ \ 

V t 2 e iy + r 2 e ix rt(e ix - e iy ) ) ' ^ ' 

The trace of this matrix is not constant, its determinant is —1 hence the spectral 
bands are not flat, thus the spectrum of the restriction of U 2 is purely absolutely 
continuous. By an analogous argument this also holds for the restriction to P ± . 
■ 

Remark that a more general periodic distribution of phases leads to ma- 
trix valued translation operators with periodic coefficients thus to non-trivial 
Hofstadter like problems. 



3 Restriction to a cylinder, transfer matrices 

Let M G N. Use the notation Z 2 M ■= Z/(2MZ) for the discrete circle of 
perimeter 2M . Consider the restriction of the model to the cylinder Z x 1*2M'- 

U(p) : l 2 (Z x Z 2M ) -> l 2 (Z x Z 2M ) 

defined by its matrix elements with respect to the canonical basis 

,v :— U(^ ljfl2 mo( i 2M);(ui,i/2 mod 2M) ■ (4) 

Remark that U(p) has the same spectral properties for some extreme cases 
as U(p), the model on the full lattice: 

Proposition 3.1 Let rt = 0. Then, for any p E CI, the spectrum ofU(p) is 
pure point. 

Proof: Similar to the proof of Proposition 2.1. I 

Proposition 3.2 Let rt / 0. Then the spectrum ofU(. ..,1,1,1,...) is purely 
absolutely continuous. 

Proof: In the proof of proposition 2.2 note that V now acts from Z 2 (Z x Z 2 m) 

to Z 2 (Z x Zm) (8) C 4 and replace T by the Fourier transform from £ 2 (Z x Zm) 
to L 2 (T x Z M ) defined by 

Tip(x,K)= ^2 ipj,ke ixj e i2 M Rk 
jeZ,JfceZ M 
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which diagonalizes the translations. Then, setting y = j^K (k G %m), the 
matrix valued multiplication operator obtained in (3) is understood as a family 
of matrix valued operators indexed over Zm- The spectral bands are not flat 
by the same argument. I 

From now on we restrict the discussion to the case 



rt / 0. 

In the following z denotes a complex number; also, unless otherwise stated, 
all indices in the second variable are to be understood mod 2M, e.g.: 



"02^2/0+1 — 1p2j,2k+l mod 2M — ^2j,2k+l[2M] ■ 

A standard approach to the spectral problem of U is the transfer matrix method. 
Though this is well known, we wish to recall the construction explicitly for the 
model at hand. 



Proposition 3.3 For z 7^ 0, q = (qi,q2,qs) £ T 3 define 
T eo {z,q) := 



T oe (z,i 



qiq2 \ 1 / z- 1 -r \ ( q 3 

q3 J t \ -r z J \ qlq 2 

(& \ 1 / z -t \ ( qjq 2 

qiq 2 J r\ t -z' 1 J\ qj 



Then 

1. Forip:Zx Z 2M it holds: 



^ U^tp v = ztp^ VfiGZx Z 2M 



^ 2j+1 ' 2k \ = T eo (z, Pe (2j,2k)) ( ^ 2j ' 2k 

1p2j+l,2k+l J e ° V ^2j,2fc+l 

and 

^' +2 ' 2fc+1 ) = T oe (z, Po (2j + l,2k + 1)) ( ^i+W 

W2j+2,2k+2 J \ W2j+l,2k+2 

2. For zeT, it holds that T oe , T eo G U (1, 1), the Lorentz group defined as 
a subset of the complex 2x2 matrices by 



U(l, 1):= Be M 2j 2(C); B*JB = J, J :- 



1 

-1 
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Proof: 

By definition of U we have for the "even" nodes: 

WW* ^ = S {pe{2j: 2k)) ( ^,2k \ = J thj+W 



{Ui>)2j,2k+1 J ' V ^2j+l,2fc+l / V ^2j,2fe+l 

and, for the "odd" nodes: 

S(p (2j + l,2fc + l))( ^i+2,2fc+2 
V ^2j+l,2fe+2 / V V2j+l,2k+l 



For a matrix 



it holds: 



5 1 = ( C- 12 ) with ^22^21 / 

<->21 <->22 



a ^ s( x )^( a )=s( x )^-( x )=s f b 



b I V v J \ y J \ b J \ a J \ y 



with 

5*22 V — ^21 1 / ' 52i V Su ~ det S 



Now 



t -r \ _ 1 / 1 -r V f t -r \ _ 1 / 1 -t 
r t J t\ -r 1 /' I r i ) ~ r It -1 



a \ /^2 \ / t —r\fq3 \ / x 



so 

fry 91 V 92 7 V r 1 J V *3 J V y 

a \ = / gig 2 \ 1 / z^ 1 -r \ ( q 3 \ / x 
V J ~\ q 3 Jt\-r z J V 7EQ2 J V b 

X \ = ( 93 V f z U M2 W & 

aj V o gig 2 J r V t -z- 1 J v o © y V y 

from which the first claim follows. 

Denote by I the identity matrix in C 2 . S is a unitary matrix if and only 

if the pullback of the quadratic form in C 4 associated with Q = ^ ^ ^ 
(blanks stand for entries) to the graph of S: { (u, Su) € C 4 , u € C 2 } is zero. 
The mapping from (x, y, a, b) to (x, 6, a, y) transforms Q to ^ ^ j V The 

pullback of the corresponding form to the graph of T eo being zero, it follows 
that T eo and, by the analogous argument, T oe , belong to the Lorentz group. 
■ 

For later use we fix the following notation 



7 



Definition 3.4 Denote by J the 2M x 2M block diagonal matrix consisting of 
M non-zero diagonal blocks equal to J and by 

U M {1, 1) := {B G M 2Mi2M (C); B*W = J}. 

the unitary group of the hermitian form defined by$. 

Note that E/jy (1, 1) is isomorphic to the classical unitary group U(M,M) of 
the hermitian form \zi\ 2 + . . . + \zm\ 2 — \zm+i\ 2 ■ ■ ■ — \z2m\ 2 ■ 

4 Relevant phases 

Because of the uniform distribution it is possible to reduce the number of rel- 
evant phases in the model to two phases per node. Before proceeding we do 
this reduction. We shall repeatedly make use of 

Lemma 4.1 Let cpi, • • • , <~p n be independent and uniformly distributed random 
variables on R/Z and let A G M m , n (Z). Then, Q x , ■ ■ ■ ,0 m defined by = Aip 
are independent and uniformly distributed if and only if Rank A is maximal. 

Proof. For k G Z m it holds 

E(e^) = E(e*^' A ^) = E{e^ A *^) = 6 At%Q . 

Thus the 9 are independent and uniformly distributed if and only if E(e*(*'^) = 
Sj: if and only if Keryl* = {0}, equivalently, if and only if Rank A is maximal. 

■ 

Proposition 4.2 There exists g : — >• T such that for p G ft the evolution 
U(p) defined by (1) is unitarily equivalent to 

D(g(p))8 on 1 2 {Z 2 ) 

where D(q) is diagonal, D(q)(j k y^ k \ = q^ k , andE> = Z7(. . . , 1, 1, 1, . . .). More- 
over, the image measure of (g>^ 2 z) 2 d I by g is ®j2dl. 

Proof: By (1), U(p) is of the form U(p) = (p)SD® (p) where flW(p) 
are diagonal, and defined by their diagonal elements: 

DW(p) 2j+h2k = p lP2 (2j, 2k), DW(p) 2j , 2k+l = Pl p 2 (2j, 2k), 

(2j + 1, 2k + 1), DW(p) 2j+h2k+1 = p 4 p 5 (2j + l,2k + 1), 

D^( P ) 2j ,2k = Ps(2j, 2k), D^{p) 2j+l ^ 2k+l = p 3 (2j, 2k), 

D^\ P ) 2j +2,2k + l = Pz{2j + 1, 2k + 1), £> (2) (p)2j+i,2fc+2 = Pe{2j + l,2fe + l). 
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Hence, U(p) is unitarily equivalent to D^> (p) DW(p)S which has the asserted 
shape. Define q = g(p) by 



q(2j + 1, 2k) := p 6 (2j + 1, 2k - l)pu> 2 (2j ) 2fc), 

g(2j, 2fc + 1) := p 6 (2j - 1, 2fc + l)pipa(2i, 2fc), 

q(2j + 2,2k + 2) := p 3 (2j + 2, 2A: + 2)p iP5 (2j + 1, 2fc + 1), 

g(2j + 1,2k + 1) := p 3 (2j, 2k)p 4 p 5 (2j + l,2k + l). 

Now, an application of Lemma 4.1 shows the q's are i.i.d. and uniformly dis- 
tributed. I 

Remark 4.3 Note that the unitary transformation just constructed is diagonal 
and thus does not affect the localization properties of the model. 

In the following, we abuse notations and call for q G T^ 2 the matrix operator 
D(q)S again U(q); same abuse for the restriction to the cylinder. 



5 Characteristic exponents 

We now define and analyze the transfer matrices and in particular the localiza- 
tion length. Consider 

U(p) = D(p)§ on l 2 (Z x Z 2M ) (5) 

with identically distributed uniformly distributed phases in T ZxZ2M , <8>zxZ 2 m c ^ 
and the cylinder set algebra. 

We use the unitary equivalence 

i 2 (ZxZ 2M )^l 2 (Zx{0,...,2M-l}) -»• 1 2 (Z;1 2 (Z 2 m)) = 1 2 (%,C 2M ) 

i/j i->- * (6) 

*j := (lpj,0,...,ll>j,2M-l)- 

Note that with the reduced phases the building blocks of the transfer ma- 
trices read with phases p,q 

f p \ 1 / -rWl \ 

/ 1 \ 1 / z -i Wg 0\ 

V o p y r V * j v o i J • 

As we shall explain below, the previous analysis leads us to deal with the 
following random dynamical system: 
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Consider the probability space defined by Q = (T 4M ) P = ® z d 4M l, and 
J 7 the cylinder set algebra. The shift 

Q-M^n, Qp(2m) :=p(2(m + 1)) (m € Z) 

is measure preserving and ergodic. For p G f2 define the following elements of 

p r = (l,pi, l,p 3 , . . . , l,P2Af-l) 
= (PO, 1,P2, 1, ■ ■ ■ ,P2M-2, 1) 
Pm = (P2M,P2M+1, • • • ,P4M-l)- 



Denote for g G T 2M the unitary diagonal matrix 



9i 



D(q) :- 



\ 



Q2M / 



(where valued matrix entries are represented by blanks) and for z / the 
2M x 2M matrices 



/ z 1 — r 



Mi(z) := 



\ 



-r z 



z 1 — r 
—rzj 



( -z- 1 



M 2 (z) :-- 



t -z' 



t -z- 



\ -t 



Define for a fixed z ^ 



A z -.n -> c/ M (i,i) 

A z (p) := D( Pl )M 2 (z)D( Pm )M l (z)D(p r ). 
Then A generates the cocycle <I> over the ergodic dynamical system 

(n,^,P,(G n ) n6Z ) 



(7) 
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defined by <£» 2 : Z x O ->• [7 M (1, 1) 

<J> z (n,p) := < I n = . 

[ i i-i(e fi p)...A- 1 (e- 1 p) n<o 

Oseledets theorem holds for <£, see [A], Theorem 3.4.11 and Remark 3.4.10 (ii): 

Definition 5.1 Let z / 0. There exists an invariant subset of full measure of 
pgfi such that the limits 

lim (<S>* z (n,p)<S> z (n,p)) 1/2n = lim (<S>* z (n,p)$ z (n,p)) 1/2H =: ^Jp) 

n— >oo n— > —dc 

exist. Denote by 7fc(p, z) A; G {1,...,2M} t/je eigenvalues of^ z {p) ar- 
ranged in decreasing order. Due to ergodicity there exists 7fc(z) > sue/? 
t/7at 7/c(p, 2) = 7&(z) on an invariant subset of full measure. The characteristic 
exponents are defined by Xk(z) := log 7^(2;). 

Due to the Lorentz symmetry of the transfer matrices for z G T we have 

Proposition 5.2 1. Let B G ?7m(1, 1)- Then for the singular values SV(B) 
it holds 

7 G 51/(5) ^ - G SF(B). 

7 

2. For Aj := log7j, 7j G 5'V(S) arranged in decreasing order it holds: 

X j+M = -Xm-j+1 Vj G {0, . . . , M}. 

Proof: We have B*$B=S. In particular det B / 0, so 7 / and 

J- 1 ^* = B^JT 1 as well as BS = IB*' 1 . Now 

det (B*B - z 2 ) = det (j-^BJ - z 2 ) = 

det ((^S)- 1 - z 2 ) = z 4M det( J B*S)- 1 det - B*B^j = 0. 

From which the two claims follow. H 
Thus we restrict our discussions to the first M non-negative Lyapunov ex- 
ponents 

Ai > A 2 > . . . > A M > 

which we shall call for simplicity "the" Lyapunov exponents in the sequel. 

We show that due to the translation invariance of the uniform distribution, 
the exponents are independent of z: 

Lemma 5.3 For any w G T, 

A wz (p) = A z (wQp), 
where w p is defined by w p2j '■= w~ 1 p2j, and w P2j+i '■= wP2j+i- 
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Proof: Write D{p m ) in (7) as D{p' l )D(p' r ), where 

Pr := (l>P2Af+l, l)P2Af+3) ■ • • , l,P4Af-l) 
Pi ■ = (P2M, l,P2M+2, 1) • • ■ ,P4M-2, 1)- 

Thus is the product of the block diagonal matrices D(pi)M2(z)D(p',) 

and D(p' r ,)Mi(z)D(p r ) whose blocks are 



a z (p, q) 



1 0\l/z -t\/10\_l/a! -9* 
p J r \ t — z^ 1 / V g ) r \ pt —pqz" 1 



r( . _ / 1 \ 1 / z" 1 -rUl \ _ 1 / z- 1 -qr 
a z{Pi<l) — ^ o p J t \ -r z J \ q J r\ -pr pqz 

For any «j£T, these matrices satisfy 

oLz(P>Q) = waiiw^P'W^q), a wz(P,q) = w~ l a l z {wp,wq), 

from which the result follows. B 
Therefore, for any fixed w € T, the matrices A z (w Q p) have the same 
distribution as A wz {p). As a consequence 

Corollary 5.4 All characteristic exponents \k{z) = \k are independent of z € 
T. 

Proof. A fc = E(log(7 fc (^p)))=E(log(7 fc (l,i©p))) = E (log( 7fc (l,p))). 



Definition 5.5 The localization length £m G [0, oo] is defined as 

£M ■= T • 

Am 

Remark 5.6 In the physics literature, see [KOK], £m in assumed to be finite for 
all parameters; a change of the asymptotic behavior as M — > oo is conjectured 
when the parameters of the model approach the critical point t = r. This 
conjecture is supported by a numerical finite size scaling method and is supposed 
to reflect the divergence of the localization length of the full system at the 
critical point. Thus a first step to support these heuristics is to prove finiteness 
of and to establish precise information of its behavior as a function of M. 

The announced equivalence to the propagation problem is the content of 
the following 

Proposition 5.7 Let U(p) be the ergodic family of unitary operators defined 
in (5) over the probability space T := T ZxZ2M , 0z,xZ 2M dl and the cylinder set 
algebra. Let f : F — > Q be defined for j € Z by 

/(P)(2j) := (P2j,0,P2j+2,l,P2j,2,P2j+2,3 ■ ■ ■ P2j,2M-2 , P2j+2,2M-1 
P2j+l,0,P2j+l,l ■ ■ ■P2j+l,2M~l) ■ 
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The image measure by f is the measure on Q and it holds 

U{p)1> = zij>^ V 2N = $*(JV, /(p))* 

for ^ defined in (6). 

Proof. The construction of / follows from Proposition (3.3). The image 
measure follows from lemma (4.1). I 

6 Finiteness of the localization length 

Using the methods exposed in [BL], see also [GM], we prove that all Lyapunov 
exponents are distinct and in particular that the localization length for the 
cylinder in finite. 

Theorem 6.1 For rt ^ 0, z G T it holds 

Ai > A 2 > • • . > Am > 0. 

Proof. We follow the strategy exposed in [BL] and prove the theorem in 
several steps making use of lemmata to be proven below. Denote by 

G := the smallest subgroup of C/m(1, 1) generated by {A(p),p G SI}. 

By lemma 6.2 

G = C/m(1,1). 

In particular it is then known : 

G is connected. 

Furthermore, see also [RS], G is isomorphic to the complex symplectic group. 
Indeed : denote by 2 the 2M x 2M block diagonal matrix consisting of M non 

zero blocks a = ( !j* ^ j ; we write: 1 = ©f f c for short; denote by 

Sp(M, C) := {B € M 2M ,2m(C); B*lB = 2} 
the complex symplectic group. From 

it follows defining C := 0f ^ ( } "jM that 

G = C/ M (l,l) = CS P (M,C)C*. 
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In order to freely use results in [BL] we shall do our argument for real matrices. 
To this end we separate real and imaginary parts and consider 

r : M 2M ,2A/(C) -> M 4 a/,4m(K) 
x = a + ib i— > 

It holds: t(x + y) = t(x) + r(y); r(xy) = r(x)r(y); t(x*) = r(x) t ; ker{r} = 
{0};detr(x) = |detx| 2 , thus 

t(C*GC) C Sp(2M,R) 
with the real symplectic group 

Sp(2M,R) := {B £ M 4 m,4m(R); B^B = 2} 

for J = ®\ M a. 

As detr(x) = | detx| 2 implies that r(x) shares its eigenvalues with x with 
the degeneracies doubled. So the Lyapunov exponents 7 defined by the r — C 
transformed products of transfer matrices are 

71 = 72 = Al > . . . > 72 P -1 = 72p = A p > . . . > 72A/-1 = 72M = Am- 

As t(C*GC) is connected one can infer from [BL] Theorem 3.4 and Exercice 
2.9 forpe {1,...,M}: 

t(C*GC) L 2p irreducible ^ 

and > =^ 72 P = X p > \ p+ \ = 72 P +i 

t(C*GC) 2p contracting J 

in particular for p = M: X]\j > 0. Now by lemma 6.3 and lemma 6.4 the group 
t(C*GC) is 2p irreducible and 2p contracting for all p G {1,...,M} so all 
Lyapunov exponents are distinct and Am > 0. I 

The following lemmata complete the proof of theorem 6.1, we use the no- 
tations introduced in the above proof. 

Lemma 6.2 

G = U M (1,1). 

Proof. By definition G C U M (1, 1) is a closed subgroup of Gl(2M,C) thus 
G is a Liegroup. By connectedness of Um{1, 1) it is sufficient to show that the 
Lie algebras q and um(T 1) coincide. Now 

u M (l, l) = {Ae M 2M>2M (C);A jk = -A k ,(-l) k ^} 

whose dimension as a real vector space equals AM 2 . 
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Denote by Dj(t) = diag(l, 1, . . . , 1, e lt , 1, . . . , 1) the unitary matrix where 
the phase sits at the j'th slot, for j = 1,2,..., 2M and use the Mj as defined 
in section 5. For zeT the matrices 

iM 2 (z)\j)(j\M 2 (z)-\ iM 1 {z)- l \j){j\M 1 {z) (8) 

belong to g, for j = 1,2,..., 2M as they are the generators of the curves 

Dj(t), M 2 {z)D j (t)M 2 1 (z), M{ 1 (z)D j (t)M 1 (z) which lie in G as 

Dj(t) = D j (t)M 2 (z)M 1 (z)(M 2 (z)M 1 (z))-\ 
M 2 (z)D j (t)M^ 1 (z) = M 2 {z)D j {t)M 1 (z){M 2 {z)M 1 (z))-\ 
M{ 1 {z)D j (t)M 1 {z) = (M 2 (z)M 1 (z))- 1 M 2 (z) J D j (t)M 1 (z). 

The generators in (8) have the same block structure as the Mj. We compute 
the relevant blocks. For iM 2 (z)\j) (j\M 2 (z)~ 1 we get 

J_ fz -t\ fl 0\fz -t\ _ i_fl -tz 
r 2 \t -z) \0 OJ \t -z) ~ r 2 [tz -t 2 

i (z -t\ /0 0\ fz -t\ i f-t 2 tz 



<- 2 \ t —zl \0 lj V t —z J r 2 \ —tz 1 

Similarly, for iMi(zy l \j) (j\M\(z) , the blocks take the form 

i f z r\ fl 0\ f z —r\ i f 1 — rz 

j? \r z) \0 OJ \-r z) ~ t 2 [rz -r 2 

i fz r\ f 0\ f z —r\ i f —r 2 rz 

1?\r z) \0 lj \-r z J *2 \-rz 1 

Now using these matrices for z ^ R and the diagonal matrix , j = 1,2 

one gets by taking suitable real linear combinations of the matrices above that, 
in both cases, the relevant blocks are generated by 

i o\ . fo o\ . / o A /o i 



0/ '" V0 II ,% V-l 01 ' VI 



For real z use the curves Dj(t), D e M 2 Dj(t)M 2 ~ 1 D- 1 ,D- 1 M{ l D j (t)M 1 D , 
with D e = © f ) and D a = © f ^ ) for u; E T which amounts to 



v W J V 1 , 

perform the change z h-» 

Taking into account the shift in the blocks and the period 2M of the indices 
in the matrices, we get that the restrictions of g and Um(1, 1) to their tridiagonal 
elements, mod 2M coincide. 

To go off the diagonals we use commutators, i.e. we exploit that X, Y G g 
implies [X, Y] G g. 
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Let A k = \k + l)(k\ + \k)(k + 1| G g, for k G Z 2 m- Considering [A,, A fc ] for 
all values of j, k, we generate a basis of all anti self-adjoint matrices that have 
non zero real matrix elements at distance two away from the diagonal (and in the 
corners, by periodicity). By commuting Ak with Aj = — G 

q, we get a basis of self-adjoint matrices with non zero purely imaginary elements 
on the same upper and lower diagonals (plus corners) only. These matrices 
correspond to the restriction of all matrices in um (1, 1) to these diagonals. 

We generalize the argument as follows: Assume we already generated a basis 
of all matrices A G um(1, 1) such that Ajk = if \ j — k\ > m, m fixed. Again, 
periodicity is implicit here. 

Let U M (1, 1) 3 Bf(m) = \ j + m)(j\ ± \ + m\. We compute 

[A j+m ,Bf(m)]=Bj(m + l). 

This way we generate all matrices A G Uw (1, 1) such Ajk = if \j — k\ > 
m + 1. 

Hence by induction, we see that q = um (1, 1), so that G = L/m(1, 1). 



Lemma 6.3 r (C*GC) = r (Sp(M, C)) is L 2p irreducible for p G {1, . . . , M}. 

Proof. Denote e i; i G {1, ... ,4M} t/je canonical basis vectors ofR 4M . 
By definition (see [BL] with adaptation to our symplectic form) 

L q := span {v G A g M 4M ; u = Met A Me 3 A ... A Me 2? -i, M G 5p(2M, M)} 

for q < 2M. Remark that the set of directions in L q corresponds to the set of 
isotropic subspaces of M 4M . 

t{C*GC) is L g -irreducible if there is no proper linear subspace V C L q 
invariant under A 9 r (C*GC). 

Consider M real numbers 

01 > 02 > . . . OjVf > 1. 

The 4M x AM diagonal matrix 

A ( 1 1 111 1 

A = diag ax, —,a 2 , ,a M , ,Oi, —,a 2 , — , • • ■ ,a M , 

V «i a 2 a M ax a 2 a M 

belongs to r (C*GC) and ex Ae3 A. . . Ae 2 g-i is an eigenvector of (A q A) n for all 
n with simple dominant eigenvalue > 1. Thus for an invariant subspace V of L q 
either eiAe 3 A. . .Ae 2(J _i G V which implies A q M(ex Ae 3 A. . .Ae 2(? _i) G V,VM, 
thus V = L q , or ei A e3 A . . . A e 2g _i G V 1 - which implies for all w G F 

= (A q M t w, ex A e 3 A . . . A e 2q -x) = (w, A q Mex A e 3 A ... A e 2q -x) 

thus y- 1 = Lq <J=^- V = {0}. Thus we conclude the claimed irreducibility for 
q = 2p, 
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Lemma 6.4 r (C*GC) = r {Sp{M, C)) is2p contracting forp G {1, . . . , 2M- 
1}- 

Proof. For any a E R \ there exist x } y G E with x 2 — y 2 = 1 such 

(x v \ 
x )< which belongs to E7 (1, 1), has eigenvalues a,l/a. Taking such 

matrices as blocks one sees that there exists an element of ?7m(1, 1) whose 
singular values are distinct: a\ > a 2 > . . ■ > clm > 1 > 1/om • • • and thus an 
element of r (C*GC) with 2M distinct singular values 

h = b 2 = oi > . . . > b 2p -i = b 2p = a p > ... > b 2 M-\ = b 2 M = a M > 0. 

Thus &2p+i/fr2p — >n^oo and it follows from proposition 2.1, p. 81 of [BL] 
that r (C*GC) is 2p contracting. ■ 

Remark 6.5 7b summarize: we have proved that if the transfer matrices gen- 
erate the complex symplectic group Sp(M,C) then the results of [BL] apply, 
i.e.: the Lyapunov spectrum is simple. The results in [BL] are stated for real 
groups only. While it is remarked in their introduction that these results should 
hold in the complex case, this seems not to be obvious to specialists in the field. 

7 Thouless formula and the mean Lyapunov 
exponent 

In this section we shall prove the announced identity in a series of lemmata. 
The quasienergy will be called z from now on. 

Theorem 7.1 Let M G N. For the first M Lyapunov exponents associated 
with U defined in Definition (5.1) it holds: 

\ Ai = - log — > - log 2 
M ^ 2 15 rt ~ 2 B 

Proof. Let z G T. Denoting by P 2 l{z) the propagator 

p 2L (z)-.n ->■ u M (i,i) 

P 2L (z)( P ) := $ z (L,p) (^(--M)" 1 
we have for m G {1, . . . , M} 

V A, = lim — log || A m P 2L (z)(p)\\ p a.e. (9) 

i 

where A m denotes the m-th exterior product; (c.f. [A], ch.3). 
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We analyse the above limit in proposition 7.2 below and show: 

M 



1 f 1 ] 

— X i = 2 / lo S \ z ~ x\dl(x) + - log - 
i Jt t 

The assertion follows by an explicit calculation proving that 

log \z — x\dl(x) = 0. 



Proposition 7.2 (Thouless formula) Let MeN,zeC\0 then 
M 

M 



1 M f 11 



Proof: Will be done in Appendix 1. I 

Remark 7.3 Remark that we prove in particular that the density of states is 
the Lebesgue measure, see lemma 9.3 below. 

7.1 Bounds on the localization length 

We now use the Thouless formula and an M independent bound on the largest 
Lyapunov exponent to derive a bound on the localization length. We remark 
that this bound is very crude and that more involved techniques should be 
established to get more detailed information; c.f. [RS] and references therein. 

First observe that a lower bound on the mean Lyapunov exponent together 
with a tight upper bound on the largest, implies a lower bound on all. 

Lemma 7.4 Let n > 0, 5 > such that VM G N, z G T 

1 M 

— > A,- > k, and Ai < k + 5, 

3=1 

then, for all j = 0, 1, . . . , M - 1, 

Xj+lZK-j^-j. (10) 

Proof: First note that \\ > jj Sf=i Aj- Thus Ai > k, which corresponds 
to (10) for j = 0. Similarly, using also the upper bound on Ai, we have for any 

1 < j < M - 1, 

j M M 

Mk<{Y j + Yj ) X k < j{* + &)+ Yl Xk 

k=l k=j+l k=j+l 
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so that 

1 M jS 

Ao 4-1 > > Xu > K . 

~~ M — j ^ k ~ M-j 



Remark: 

In view of localization properties, the estimate is useful only if 

K>(M-l)d. (11) 

We now estimate the cocycle to derive an upper bound on the largest Lya- 
punov exponent, which is uniform in the quasienergy and width of the strip M. 



Proposition 7.5 Let M e N 

1. For the generator of the cocycle defined in (7) it holds 

||A(p)||<I(l + r)(l + t); 

2. it follows: 2Ai < log (i) + log ((1 + r)(l + *)). 

3. 7/7ere ex/'sts a c > suc/7 tftat for M € N /'t holds: 



dist(r, {0, 1}) < e 



-cAf 



Cm = v 1 - < 



A M " Iog(i)-(M-l)log((l + r)(l + t))" 

Proof. The estimate on A follows from its definition. The estimate on Ai is 
obtained using the equality (9). Finally, from the estimate (10) it follows 

1 2 

£m = i— < 



Am " log (i) — (M — 1) log ((1 + r)(l + *)) 
id 

from the critical point ^= because of the singularity of logl/rt. 



The bound is symmetric around t = r = -4= and finite for r sufficiently away 

V 2 



8 Spectral Localization 

We follow the strategy which was successfully employed for the case of one di- 
mensional Schrodinger operators: polynomial boundedness of generalized eigen- 
functions, positivity of the Lyapunov exponent and spectral averaging. We lean 
on the work of [BHJ, HJS]. Our result is: 



19 



Theorem 8.1 Let M G N, rt / 0. Then, the Chalker Coddington model 
on the cylinder exhibits spectral localization throughout the spectrum, almost 
surely. More precisely, 

1. the almost sure : spectrum E, continuous spectrum E c and pure point 
spectrum E pp ofU(p) satisfy 

E = E pp = T and E c = 0; 

2. the eigenfunctions decay exponentially, almost surely. 
Proof: We prove the theorem in Appendix 2. 

9 Appendix 1 

We follow the strategy of [CS2] and first prove the lower bound 

1 M f 11 

M Yl X ^ z ) > 2 / log |z - x\dl{x) + - log - - log \z\ (12) 

for / z G C \ T which follows from Lemma 9.3 equation (15) below in the 
limit L — > oo. 

Lemma 9.1 Denote U D the unitary defined by restriction ofU to 
l 2 ({— 2L, . . . , 2L}, 1 2 (Ij2m)) with reflecting boundary conditions: the scatter- 
ing picture for the links which are incoming to walls at — (2L + 1) and 2L + 1 
reads 

U°e-2L,2k+l = e-2L,2fc+2, U°e2L,2k = e-2L,2k+l- 

For z eC let 

F z := {if; £ l 2 (%2M);ip2k+i = zip 2 k+2, k G Z M } , 

G z := {if; G 1 2 (Z 2 m); zip2k+i = ip2k, k G Z M } 
and denote Qp the orthogonal projection to a subspace F. It holds: 

z is an eigenvalue ofU D 

$>2L = P 2 l(z)^-2L and f_ 2 L G F z and G G z 
^Ker (Q Gi P2L{z)Q F ) ^{0} 

Proof. It holds 

U D if) = zif)^ 

lfj-2L,2k+l = Zlp-2L,2k+2 and Zlp2L,2k+l = 1p2L,2k 
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so 

*-2L G and G G 2 . 

The identity 

holds by construction of the transfer matrices so 

U D i> = zip^ Q G xP 2L (z)Q Fz t> = 0. 



Lemma 9.2 Denote the "even" subspace ofl 2 (7L 2 M) by 

E := span{e 2k ; k G %m}. 

For z / there exist invertible operators V z , W z on 1 2 (Z 2 m) such that W Z (E) = 
F z and V Z (E) = Gj such that 

1. 

z is an eigenvalue ofU D 

det {QeV z ~ 1 P 2 l(z)W z Q e ) = 

where we understand the determinant to apply to the restriction to E. 
2. For z 7^ 0; {zi, . . . , z^ l+ ^ 2 m} the eigenvalues ofU D it holds: 

\z\^ M \ det (QeV~ 1 P 2 lW z Qe) I = j^ri£ +1)2M \z - *|. 

(13) 

Proof. Fix 

In the following Nj, Dj denote generic, z independent matrices whose precise 
values may change from line to line. The Dj are diagonal. 
The transfer matrix A z defined in (7) is of the form 

A z = - (z 2 D 1 Q + z~ 2 D 2 Q E + zN x + N 2 + z' 1 ^) 
rt ' 

where O denotes the "odd" subspace defined by O + E = I 2 (Ij 2 m)- Thus 

4L-1 

(rt) 2L P 2L = z 4L DxQ + z- 4L D 2 Q E + ^ z j Nj. 

j=-4L+l 

Note that 

F z = span | ~~ (ze 2k+1 + e 2k+2 ) ;k € Z M 



21 



G z = span | (e 2 fc + z 1 e 2 fc+i) ; k G Z M | • 
On 1 2 (Z 2 m) define the operators 

:= ^ ke 2fc+ i + e 2 fe + 2)(e 2 fc+2| + I - e 2 fc+i + ^~ 1 e 2 fc+2)(e2fc+il 

:= ^= ^ | - e 2fc+ i + ze 2k )(e 2k \ + \e 2k + z~ 1 e 2k+1 )(e 2 k+i\- 
Then W z Qe = Qf z W z and 14Q.E = Qq±V z . Moreover one checks that 

w 2 = i, V7 1 = in4# (14) 

with K := Y.k€L M |e2fc+i)(e2fc| + |e 2 fc) <e 2 fc+i| . It follows: 

z is eigenvalue of U D ker {QeV z 1 P 2L W z Q e ) / {0}. 

Now 

QeV z X P2lW z Q e = 
^ ^2 l e2fc ) \ ~ e 2fc+i + = e 2fc, -P2L (^e 2m +i + e 2m+2 )\ (e 2m+2 \ = 

k,m 

-±—l z ^D 1 Qo + z~ 4L ~ 1 D 2 Q E + V 
Multiplication by z 4L+1 implies that for some aj G C 

(8L+2)M 

z^ M det(Q E V z ^P 2L W z Q E ) = ^ LM £ *V 

Z?i is unitary thus, in particular, \ci( 8 l +2 )m\ = 1- z( 4L+1 ) M det . . . being a poly- 
nomial of degree (8L + 2)M whose leading coefficient has modulus (rt)~ 2LM 
and which is zero on the (4L + 1)2M eigenvalues of C/ D the formula for the 
determinant follows. 

■ 

We now prove convergence of the finite volume (L < oo) density of states 
^ as L — > oo to a non random measure: the density of state. Then we show 
that this measure equals the Lebesgue measure. 

Lemma 9.3 Denote the measure defined by 

(AL + \)2M trHUD) =; </ 6 C < T »' 

77)en 
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1. 



M vaguely 



the Lebesgue measure on T. 
For M G N there exists c M > such that for all L G N, / z G C \ T 



~log||A-P 2i || 

> ~log~ + (2 + ~)|log|z-*|d//f(x) (15) 

Proof. 1.: We first prove the existence of a nonrandom limit measure. The 
first step consists in showing that p a.e, 



lim 

L— >co 



[ f(U D {p))d^ = - {E «eo |0 , f(U)eo,o)) + E « ei)1) f(U)e 1;1 )) 
Jt 4 



E((ei,o,/(C/)ei,o))+E((eo,i,/(C/)e ,i»}=: / fdn 



for all / G C(T). This follows from a classical argument based on ergodicity, 
separability of C(T) and that U - U D has norm and rank uniformly bounded 
in L, see e.g. [Jl] for the details for the unitary case. 

In order to identify recall that the normalized Lebesgue measure dl on 
T is uniquely characterized by : 

p n dl = 5 n fi n G Z 

T 

where p n (x) := x n . Consider the space of loops of euclidean length n starting 

at (0,0) : 

r (0 ,o) = {7 = {0, • • • , n} -> {-2L, . . . , 2L} x Z 2M , 7 (0) = 7 (n) = (0, 0)} . 
Then because of the structure of U 

{eo,o,Pn(U)e 0t o) = ^ (e ,o, Ue^) . . . (e 7 ( n _ 1 ), Ue ,o). 
7er (0i0) 

Now (e 7 (j), C/(p)e 7 (j +1 )) = l(p)t a r^ for some a,/3 > and / a uniformly 
distributed random variable. Thus E ((eo,o,Pn(^) e o,o)) = ^n,0- Applying the 
same argument to (e^i, f(U)e\ : i), (eo,i ... we conclude: 

d/^ = dZ 
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2.: By formula (13): 
1 



4LM 



logjdet (QeV- 1 P 2L W z Q e ) 



1 , 1 

2 



z\ < 



log i + (2 + ^) jf log |z - x|d/if (x) - M + ^ log 

^ log || a a/ (QeVt^l^Qb) || < 

~ log || A M P 2i || + 1 (l og || /\ M Q eV -1\\ + log || A M W z Qe\ 



—•cm 

where we used the identity 

det QeAQe = <e A . . . A e 2M - 2 , A M Ae A ... A e 2M - 2 } • 
From this the claim follows: 



We turn now to the proof of the opposite inequality: 

M 

M 



1 f 11 

^ A,(z) < 2 / log \z - x\dl{x) + - log - - log \z\ (16) 

for ^ z G C \ T : 



Proposition 9.4 Suppose that for any choice of sets of vectors |d ,d 2 , ■ ■ ■ , ^ 2 m- 2 } 
and {c^q", dj, . . . , d 2 f A/- 2 } t ^ ie "°dd" subspace O 

lim sup — — ^ — x 

M(4L + 1) 

x log I ((e + d ) A . . . A (e 2M -2 + d^ M _ 2 ), 

/\ M (V z - 1 P 2L (z)W z )(e + 4) A ... A (e 2M - 2 + 4m- 2 )> I 

11 f 

<-ln h 2 / log b - z\ dl(x) - log \z\ (17) 

2 rt J T 

then, for all ^ z G C \ T, 

M 

M 



1 M 11/" 

— : V Aj(z) < - log — + 2 / log \x - z\ dl(x) - log \z\ 
M t—' 2 rt J j 



Proof: The vectors of the form {(e +d )f\. . -A(e 2 M- 2 +d 2 A/- 2 ); d , . . .,d 2 M-2 G 
0} span A M C 2M . On the other hand, given any spanning sets Si and S 2 in 
A M C 2M , the mapping || • || 5 defined by 

s= sup 1(0,^)1 
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defines a norm over the algebra of operators in AC. It follows that there 
exists c > 0, which depends on Si and S 2 , such that for any matrix A, \\A\\s > 
c[|A||, hence that 

1 M 1 1 r 

V Xi(z) < - log - + 2 / log |C - z\ dl{() - log \z\ . 

■ 

We now prove that the inequality (17) is satisfied. This will be achieved in 
two steps. In order to keep track of the L dependence denote by the former 
U D . Now reinterpret the left hand side of (17) as the characteristic polynomial 
of a deformation of denoted V^; more precisely: we aim at equation (18) 
below. The problem is then reduced to the proof of the weak convergence of 
the associated sequence of counting measures iy^z) towards fi . 

9.1 Deformation of C/f 

Let L in N and define the matrix Vf +1 on 1 2 ({-2L - 2, . . . , 2L + 2}, l 2 (Z 2M )) 
by: V4> € 1 2 ({-2L - 2, . . . , 2L + 2}, l 2 (Z 2M )), 

M-l 

(Vl+1^) 2 L+2,2* = Yl B 2k,2l^2L,2l 
1=0 

(VL +1 lp) 2 L+l,k = 4>2L+l,k 
M-l 

(^ / L+lV')2L,2fe+l = ^ C 2k+l,2l+l^L+2,2l+l , 
1=0 

M-l 

(VL+i1p)-2L,2k = Yl B 2k,2l^~2L-2,2l 
1=0 

(V£ +1 1p)-2L-l,k = ll>-2L-l,k 
M-l 

(V^ +1 Tp)^2L-2,2k+l = Y C< 2fc+l,2/+lV'-2L,2Z+l , 

1=0 

with the same reflecting boundary conditions as U£ +1 and 
(Vl+iVO^j" = (^F+i'P) ^ ^ or an y va l ue s of (/x, v ) which were not described pre- 
viously. The matrix V^_ 1 is a deformation of the matrix U£ +l , but its structure 
remains close to the structure of U^ +1 . Note that span{e2L+i,k, e~2L-l,j',j, k G 
%2m} belongs to ker (V^_ 1 — I). For ^ an eigenvector of V® +1 associated with 
the eigenvalue z, 

vf +1 v = 
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This implies that either tp G span{e2L+i,k, £-2L-i,j', j, k G Z 2 m} and z = 1, 
or tp-2L-2 £ and V2L+2 G and 



1p-2L 



P2L(z)lp-2L 

A + {z)^2L = z^QeB+Qe^l + zQoC+^QoihL 



2L-2 



QeB QeiP-2L-2 + zQoC Qo^-2L-2 



The transfer matrices ^4 + (^) and A~(z) are deformations of the matrices A z . 
This construction is useful to establish the following lemma. In the following, 
z will be fixed as a parameter. 

Lemma 9.5 Let (d 2 f /c )fc g {o,...,M-i} ar >d (d 2k )ke{o,...,M-^} two families of vectors 
belonging to the "odd" subspace O. These families are the columns of two 
corresponding matrices denoted D + and D~ respectively. Assume that for z 7^ 
0, max(||D _ ||, ||-D + ||) < max(|^|, \z\) and consider the matrixV® parametrized 



by z 



CI 



z + TD + 

z{\- z- l D + T)- 1 
z + z 2 D~* 



(z- 1 -D~ 



with 



1 

1 



1 \ 



1 oy 



Then, 



((eo + d ) A . . . A (e 2 M-2 + d 2M _ 2 ), 

A M (V- 1 P 2L (z)W z )(e + d+) A ... A (e 2M -2 + d+ M _ 2 )) 
= (e A . . . A e 2 M-2, A^-^-P^A+W^eo A ... A e 2M -2> • (18) 

Proof: By (14): W 2 = I, V~ l = KV Z K , z ^ 0. It follows for all G 

{0,...,M-1} : 

A + W^e 2fc = W*(e 2fc + d+) 

A-VT 1 *^ = V- u (e 2k + d- k ) . 
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Given z, D + and D and the associated matrix V® +1 , we consider the 
corresponding eigenvalue problem: 

v£> +1 i> = sty 

The complex number z' is an eigenvalue of V® +1 iff 

(z' - l) 4M (e A ... A e 2 M-2, A M (V^ 1 A~ (z')P 2L (z')At (z')W z ,)e A ... A e 2M ^) = 

where Af(z') = z'- x Q E B±Q E + /QoC^Qo- 

Once multiplied by z '( 4L + 5 ) M i t ne left-hand side is a polynomial of degree 
2M(AL + 5) in z' . Following the Thouless argument, we get for the logarithm 
of the modulus divided by AML: 

\ log A - log \z'\ + (2 + ~) f log |x - z'\ dvl\ z {x) + 0{\) , 

where the family of measures Vjf z are supported on some closed ball B(0,R z ), 
due to the fact that sup L \\Uf -V®\\ < oo. Note that if z' = z, Af(z) = A + 
and A~(z) = A~ . 

9.2 End of proof of inequality (17) 

We split the proof in the two following lemmas, whose proof is an adaptation 
of the argument given in [CS2]. 

Lemma 9.6 If {vl)l&$ ar *d p- are measures supported on B(0,R) for some 
R > 0, and if (ul) converge weakly to p, then for any z G C 

/ log|C- ^i(C) < / log|C-*|dM(C) • 

JT JT 

Proof: Given zeC, let f e be defined by: 

f / 6 (C) = Iog|C-^|if|C-2|>e 
\ / 6 (C) = Iog|€|if|C-«|<e 

Since the support is compact, 

lim / / e (C)«MC) = ( feiOMQ ■ 
On the other hand, for any C in T, 

logic - ^1 < /.(C) 
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so that: 

limsup / log|C - A dvUC) < limsup / fe{Qdv L {0 
L^too JT L^too JT 

= I fe(C)dKC). 
JT 

The result follows by monotone convergence theorem when e goes to zero. 
■ 

Remark: Let us note that the *-algebra of trigonometric polynomials F defined 
by: 

T = {feC(B(0,R));f(r,9)= £ a klM r k ^e^ 9 , N G N , fci G N , fc 2 G Z} 

fci + |fc 2 |<Af 

separates points and contains the constants. Its closure under the supremum 
norm is C(B(0, R)). The weak convergence of the measures is equivalent to 
have for all / in J 7 , 



lim / /(C)<MC) = / /(C) ^(0 

Jb(om) Jb(o,R) 



B(0,R) JB(0,R) 



Lemma 9.7 As a Borel measure on C, the sequence of measures {vff z ) para- 
metrized by z, M converges almost surely weakly to dl as L tends to infinity. 

Proof: Let (rj i2 e^.^)^^ 4i+1 ^ and (e x i )^/ 4L+1 ^ be the eigenvalues of the 
problems with reflecting boundary conditions for V and Ud which correspond 
respectively to the modified and unmodified "potentials". We have that: 



M 1 \ ~> r 

L > z 2M(4L + 1) V ' 
1 

2M(4L + 1) 



U M = 1 V x 

° ; ' ' 1 /. -)- 1) / / e 



These measures are supported on some B(0,R z ). We will drop the z subscript 
in the sequel. Since we already know that {pL 1 ^) converges almost surely weakly 
to dl, we only need to show that for any nonnegative integer k\ and any integer 

2M(4L+1) 

lim — V e ik2X > - r kl e ik2 ^ = . 

L^oo 2MUL + 1) ^ 3 
j=i 
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Actually, it is enough to prove it for nonnegative integers k\, ki- Let us fix such 
a couple (fci, fo) ar) d decompose the term on the left-hand side as follows: 



2M(4L+1) 

* V" e <teAi_ fci ifa^ = Tl (L) + T 2 (L) 

2M(4L + 1) ^ J 



2Af(4L+l) 

Where r '< £ > = 2M|4iTT) g tf-fV* 

r »< £ > = w(4bi) Tr(£/ «- Ffa) ' 

If = min(A;i, £2) and I = max(k\, A^), we have 

2M(4X+1) fe 2M(4L+1) 

\Ti(L)\ < ^ r V r fe |ri-- fc -l| < ^ , V Ir^-ll . 

1 v ;| -2M(4L + 1) ^ lj l_ 2M(4L + l) ~l J 



Following [CS2], we first prove that 

lim 

L^oo 2M(AL + 1) 



^ 2M(4L+1) 

lim > |r, — 1| = . 



We know that there exists two orthonormal bases (</>j)j^ 4L+1 ^ and 
such that: 

2A/(4L+1) 

vf= E «,-(y?)i^-x^i , 

where (fJ-jiVf)) are the singular values of the operator V^ 5 . Actually, Hj{V^) = 
Tj and we assume them to be ordered: /Xj + i(V^) > Hj{V®) > 0. Note that: 

0*?(VJ?); j E {1, . . . , 2M(4L + 1)}} = a{VE *V° ) \ {0}. 

Since for each j e {1, . . . , 2M(AL + 1)}, Hj(U£) = 1 we deduce from the 
remark following Theorem 1.20 in [S] that: 

2M(4L+1) 2M(4L+1) 2M(4L+1) 

E i»-i-ii= E m v l) -h<pE)\< E n<y?-vE)- 
j=i j=i j=i 

Since V^P — has rank and norm uniformly bounded in L, we obtain that: 

lim Ti(L) = . 

L— ¥00 
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The term T2(L) will be treated in a similar way. The operator — V® has 
rank and norm uniformly bounded in L. This implies that for all integer k%, 
Uf k2 — V® k2 has also rank and norm uniformly bounded in L. So, 

lim T 2 (L) = , 

L— »oo 

which concludes the proof. I 

The above lemmata together with equation (18) establish the inequality 
(17) which implies (16). We finish with the proof of the Thouless formula on 

T: 

Lemma 9.8 For all zeT, 

1 M 1 1 f 

— \i{z) = - log — + 2 log \x - z\ dl{x) 

Proof: We note with [CS1] that lim^oo i log || A m (P 2 l)(-z)|| is subharmonic 
in C \ {0} and J T log \z — x\dl(x) subharmonic on C\T. The two exceptional 
sets are of measure zero in C, these quantities must agree everywhere. I 

Remark 9.9 We note that the above proof does not depend on the specific 
form of the density of states. 

10 Appendix 2 

Now we prove theorem 8.1 in several steps. 

By theorem 6.1 the localization length is finite for all values of the param- 
eters. Note that the spectrum is characterized by the existence of generalized 
eigenfunctions: 

Suppose that the support of E p {-), the spectral resolution of U(p), is the 
whole circle T. 

Proposition 10.1 For M G N, p E Q, the spectrum ofU(p) is the closure of 
the set 

S p = {z G T; U{p)4> = zcf) has a non-trivial polynomially bounded solution} 
and E P (T \ S p ) = 0. 

Proof: The stated behaviour at infinity of the generalized eigenvectors 
and the spectrum of U{p) are related by Sh'nol's Theorem. This well known 
deterministic fact for self-adjoint operators was proven in [BHJ] to hold in the 
unitary setup for band matrices on Z 2 (Z). It is straightforward to check that the 
result holds for band matrices on l 2 (Z,C 2M ), with M finite. I 

Secondly we prove the existence of a finite cyclic subspace: 
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Lemma 10.2 Let M e N,rt / 0. Denote I := {0} x Z 2 M- The vectors 
{e M ; /i e Jo} span a cyclic subspace oft 2 (Z x Z2m)- 

Proof: The only non vanishing elements in U are the blocks given in equa- 
tion (1). Denoting generically the elements of S by 



S 



a f3 
7 5 



and observing that U„ I = U v ^ we have 

U(2j+l,2k);{2j,2k) U( 2 j+l,2k);(2j+l,2k+l) ] _ f a $ 
U(2j,2k+l);(2j,2k) ^(2i,2fe+l);(2j+l,2fc+l) / \ 7 ^ 

^ (2j+2,2fc+2);(2j+2,2fc+l) ^(2,j+2,2fc+2);(2j+l,2fc+2) 
U (2j+l,2fc+l);(2j+2,2fc+l) ^ (2j+l,2fc+l);(2j+l,2fc+2) 



U (2j,2k);(2j+l,2k) U (2j,2k);(2j,2k+l) ) _ / '> ' 



and 



\ ^(2j+l,2k+l);(2j+l,2k) ^(2j+l,2fc+l);(2j,2fc+l) J V P $ 

(^2j+2,2k+l;2j+2,2k+2 ^2j+2,2k+l;2j+l,2k+l J 
^2j+l,2fc+2;2j+2,2fc+2 ^2,j+l,2fc+2;2j+l,2£;+l / 

Computing Ue( 0t 2k) = ae (i,2k) + 7 e (o,2fc+i) ar| d the corresponding expressions 

for f/^ 1 e(i ) 2fc),f/e (0; 2A : +i),^ 1 e(_i i 2fc+i) we infer: 

e (l,2fc) = — (^e (0] 2fc) - 7 e (0,2fe+l)) 

P 7 TT -i 

e (l,2fc+l) — ~ e (0,2k) - = u e (0,2fc+l) 

e (-l,2fc+l) = ~ (^ e (0,2fc+l) - ae( 0i 2fc+2)) 

- 6 °TT-l 
e (-l,2fc+2) — - e (0,2fc+l) _ -= u e (0,2fc+2)- 

Thus vectors with indices in {±1} x Z2A/ belong to the subspace generated by 
L r=tl (/o)- The lemma follows by induction. I 

Let / = {0, 1} x Z 2M , H = T z4M \ 7 , P = 8) fe6Z 4M VJ dl, p = {pj}^^ G 
fi, and 0/ = {^j}je/- We shall use the notation Q, 3 p = (p,Qj). 
Denote 

AmM := ^lim (^log(|| A M P 2 L (*)(p)l|) - ^l°g(ll ^ ^WWII)), 
if the limit exists. 
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By construction, 5.1, it holds for almost every p 

Am = Am(p) 

By definition Xm{p, z) is independent of the finitely many 0/, if p = (p, Qj). 
By theorem 6.1 there exists Q(z) C 12 with F(Q,(z)) = 1 such that for any 

z G T\R 

Am((p, 0/), 2) = Am > 0, 

for all 9j G 0/ and all p G f2(z). We can apply Fubini to the measure P x dl to 
get the existence of S7o £ ^ with P($^o) = 1 such that for every p £ Qq there 
is Bp G T with l(Bp) = and 

\ M ((p,Qi),z) > for all % G 0/, and all z G % c . (19) 

Then we show that for p G Oo, -Bp*" 7 is a support of the spectral resolution 
of (7((p,9j)) for almost every 0j G 0/ w.r.t. d^l/ on 

For any fixed j G /, we introduce the spectral measures up associated with 
U(p) = f T x dE p (x) defined for all Borel sets A G T by 

l4(A) = {e j \E p (A)\e j ). 

Since U(p) = D(p)S, where D(p) is diagonal, the variation of a random phase 
at one site is described by a rank one perturbation. More precisely, dropping 
the variable p temporarily, we define D by taking 9j = 1 in the definition of D: 

D = D + leyXeyKl - 0,) = e l ° s(§ ^ ){e ^D, 
so that, with the obvious notations, 

U = D§ = e lo s(^j)l e J>< e J I f/. 

The unitary version of the spectral averaging formula, see [C] and [B], reads in 
our case: for any / G L 1 ^), 

[ dl(9j) [ f(xWL ) (x)= [ f(x) dl(x). 

Applied to / = xb ¥ i the characteristic function of Bp, this yields 

= l(Bp)= f M ^ ei) (%)dZ(^). (20) 

Consequently, 

A*/pQ \{Bp) = 0, for every 9% G 0/, k / j and Lebesgue-a.e. 9j, 
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Therefore, for all p G Qq, there exists Jp C T^l s.t. l{Jp C ) = and 

8/C Jp^/iJ_ e/) (%)=0, VjGl. (21) 

Now fix p E S7o and 0/ C Jp and consider p = (p, Qj). By Lemma 10.2 and 
(21) we deduce that E p (Bp) = 0. If S p is the set from Sh'nol's Theorem 10.1, 
then the set S p n Bp C is a support for E p (-). 

Now take z G S p H Bp . By Theorem 10.1, U(p)ip = zip has a non-trivial 
polynomially bounded solution if). On the other hand, by (19), \m{p,z) > 0. 
Thus, by Osceledec's Theorem, every solution which is polynomially bounded 
necessarily has to decay exponentially both at +oo and — oo, and therefore it is 
an eigenfunction of U(p). In other words, every z G S p D Bp- C is an eigenvalue 
of U(p), hence S p D Bp C is countable. Therefore E p (-) has countable support 
thus U(p) has pure point spectrum. With 

:= {(P,{0j}jei) s -t- PG^o, {%} ie /CG/ C Jp}, 

we have 

p6O a c {U{p)) = %. (22) 
Also, from l(Jp C ) = we have 

(® j6J cK)(J ? ) = (® i6/ d/)(Tl 7 l) = 1. (23) 
As P(H ) = 1, we conclude from (22) and (23) that 

F(a c (U(p)) = 0) > P(J2o) = / cflP(p)(® i6J cfl)(J ? ) = 1, 

JClo 

which proves that U(p) has almost surely pure point spectrum. The fact that 
the support of the density of state coincides with the almost sure spectrum, see 
[Jl], shows that S pp = T. 

We finally show that almost surely all eigenfunctions decay exponentially. 
Note that we actually have shown above that the event "all eigenvectors of 
U(p) decay at the rate of the smallest Lyapunov exponent" has probability one, 
since this is true for all p G firj- Measurability of this event was proven for the 
case of ergodic one-dimensional Schrodinger operators by Kotani and Simon in 
Theorem A.l of [KS]. The proof of this fact provided in [KS] carries over to 
the CC model as well. It is enough to note that, due to Lemma 10.2, we may 
use p p = YljEi A as spectral measures in their argument. I 
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